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^ ■ ABSTRACT 

, The hnearized relativistic magnetohydrodynamic (RMHD) equations describing a 

uniform axially magnetized cylindrical relativistic jet spine embedded in a uniform axi- 
I ally magnetized relativistically moving sheath are derived. The displacement current is 

retained in the equations so that effects associated with Alfven wave propagation near 
light speed can be studied. A dispersion relation for the normal modes is obtained. An- 
' alytical solutions for the normal modes in the low and high frequency limits are found 

and a general stability condition is determined. A trans-Alfvenic and even a super- 



O ' Alfenic relativistic jet spine can be stable to velocity shear driven Kelvin-Helmholtz 

^ I modes. The resonance condition for maximum growth of the normal modes is obtained 

^ ' in the kinetically and magnetically dominated regimes. Numerical solution of the dis- 

persion relation verifies the analytical solutions and is used to study the regime of high 
sound and Alfven speeds. 
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1. Introduction 



Relativistic jets are associated with active galactic nuclei and quasars (AGN), with black 
I hole binary star systems (microquasars), and are thought responsible for the gamma-ray bursts 

. (GRBs). In microquasar and AGN jets proper motions of intensity enhancements show mildly 

superluminal for the microquasar jets ~ 1.2 c (Mirabel & Rodriquez 1999), range from subluminal 
(<C c) to superluminal (< 6 c) along the M 87 jet (Biretta et al. 1995, 1999), are up to ~ 25 c along 
^ ' the 3C345 jet (Zensus et al. 1995; Steffen et al. 1995), and have inferred Lorentz factors 7 > 100 

in the GRBs (e.g., Piran 2005). The observed proper motions along microquasar and AGN jets 
imply speeds from ~ 0.9 c up to ~ 0.999 c, and the speeds inferred for the GRBs are ~ 0.99999 c. 

Jets at the larger scales may be kinetically dominated and contain relatively weak magnetic 
fields, e.g., equipartition between magnetic and gas pressure or less, but the possibility of much 
stronger magnetic fields exists close to the acceleration and collimation region. Here general rel- 
ativistic magnetohydrodynamic (GRMHD) simulations of jet formation (e.g., Koide et al. 2000; 
Nishikawa et al. 2005; De Vilhers, Hawley & Krolik 2003; De Villiers et al. 2005; Hawley & Krolik 
2006; McKinney 2006; Mizuno et al. 2006) and earlier theoretical work (e.g., Lovelace 1976; Bland- 
ford 1976; Blandford & Znajek 1977; Blandford &; Payne 1982) invoke strong magnetic fields. In 
addition to strong magnetic fields, GRMHD simulation studies of jet formation indicate that highly 
collimated high speed jets driven by the magnetic fields threading the ergosphere may themselves 
reside within a broader wind or sheath outflow driven by the magnetic fields anchored in the ac- 
cretion disk (e.g., McKinney 2006; Hawley & Krolik 2006; Mizuno et al. 2006). This configuration 
might additionally be surrounded by a less collimated accretion disk wind from the hot corona 
(e.g., Nishikawa et al. 2005). 
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That relativistic jets may have jet-wind structure is indicated by recent observations of high 
speed winds in several QSO's with speeds, ~ 0.1 — 0.4c, (Chartas, Brandt & Gahagher 2002, 
Chartas et al. 2003; Pounds et al. 2003a; Pounds et al. 2003b; Reeves, O'Brien & Ward 2003). Other 
observational evidence such as limb brightening has been interpreted as evidence for a slower external 
sheath flow surrounding a faster jet spine, e.g., Mkn 501 (Giroletti et al. 2004), M 87 (Perlman et al. 
2001), and a few other radio galaxy jets (e.g.. Swain, Bridle & Baum 1998; Giovannini et al. 2001). 
Additional circumstantial evidence such as the requirement for large Lorentz factors suggested by 
the TeV BL Lacs when contrasted with much slower observed motions suggests the presence of 
a spine-sheath morphology (Ghisellini, Tavecchio & Chiaberge 2005). At hundreds of kiloparsec 
scales Sicmignowska ct al. (2007) have proposed a two component (spine-sheath) model to explain 
the broad-band emission from the PKS 1127-145 jet. A spine-sheath jet structure has been proposed 
based on theoretical arguments (e.g., Sol et al. 1989; Henri & Pelletier 1991; Laing 1996; Meier 
2003). Similar type structure has been investigated in the context of GRB jets (e.g., Rossi, Lazzati 
k Rees 2002; Lazzatti & Begelman 2005; Zhang, Wooseley & MacFadyen 2003; Zhang, Woosley & 
Heger 2004; Morsony, Lazzati & Begelman 2006). 

In order to study the effect of strong magnetic fields and the effect of a moving wind or 
sheath around a jet or jet spine, I begin by adopting a simple system with no radial dependence of 
quantities inside the jet spine and no radial dependence of quantities outside the jet in the sheath. 
This "top hat" configuration with magnetic fields parallel to the flow can be described exactly 
by the linearized relativistic magnetohydrodynamic (RMHD) equations. This system with no 
magnetic and flow helicity is stable to current driven (CD) modes of instability (Istomin & Pariev 
1994, 1996; Lyubarskii 1999). However, this system can be unstable to Kelvin-Helmholtz (KH) 
modes of instability (Hardee 2004). This approach allows us to look at the potential KH modes 
without complications arising from coexisting CD modes (see Baty, Keppens Sz. Compte 2004) and 
predictions can be verified by numerical simulations (Mizuno, Hardee Sz Nishikawa 2006). 

This paper is organized as follows. In §2, I present the dispersion relation arising from a 
normal mode analysis of the linearized RMHD equations. Analytical approximate solutions to 
the dispersion relation for various limiting cases are given in §3. I verify the analytical solution 
through numerical solution of the dispersion relation in §4. I summarize the stability results in 
§5 and discuss the applicability of the present results in §6. Derivation of the linearized RMHD 
equations is shown in Appendix A, derivation of the normal mode dispersion relation is presented 
in Appendix B, and derivation of the analytical solutions is shown in Appendix C. 

2. The RMHD Normal Mode Dispersion Relation 

Let us analyze the stability of a spine-sheath system by modeling the jet spine as a cylinder 
of radius R, having a uniform proper density, pj, a uniform axial magnetic field, Bj = Bj^z, and a 
uniform velocity, uj = uj^z- The external sheath is assumed to have a uniform proper density, pe, a 
uniform axial magnetic field, Bg = Bg^z, a uniform velocity Ug = Ug^z, and extends to infinity. The 
sheath velocity corresponds to an outflow around the central spine if Ue,z > or represents backflow 
when Ug^z < 0. The jet spine is established in static total pressure balance with the external sheath 
where the total static uniform pressure is P* = Pg + B^/Stt = P* = Pj + Bj/Stt, and the initial 
equilibrium satisfies the zeroth order equations. Formally, the assumption of an infinite sheath 
means that a dispersion relation could be derived in the reference frame of the sheath with results 
transformed to the source/observer reference frame. However, it is not much more difficult to derive 
a dispersion relation in the source/observer frame in which analytical solutions to the dispersion 
relation take on simple revealing forms. Additionally, this approach lends itself to modeling the 
propagation and appearance of jet structures viewed in the source/observer frame, e.g., helical 
structures in the 3C 120 jet (Hardee, Walker Sz Gomez 2005). 

The general approach to analyzing the time dependent properties of this system is to linearize 
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the ideal RMHD and Maxwell equations, where the density, velocity, pressure and magnetic field 
are written as p = po + Pi > v = u + vi (we use vq = u for notational reasons) , P = Pq + Pi, and 
B = Bo + Bi, where subscript 1 refers to a perturbation to the equilibrium quantity with subscript 
0. Additionally, the Lorentz factor 7^ = (70 + 71)^ — To + ^TqU • vi/c^ where 71 = 7gU • vi/c^. 
The linearization is shown in Appendix A. In cylindrical geometry a random perturbation pi, vi 
Bi and Pi can be considered to consist of Fourier components of the form 



fi{r,(f),z,t) = fi{r) exp[i{kz ± ncp - tot)] 



(1) 



where flow is along the z-axis, and r is in the radial direction with the flow bounded by r = i?. In 

cylindrical geometry n is an integer azimuthal wavenumber, for n > waves propagate at an angle 
to the flow direction, and +n and — n give wave propagation in the clockwise and counter-clockwise 
sense, respectively, when viewed in the flow direction. In equation (1) n = 0, 1, 2, 3, 4, etc. 
correspond to pinching, helical, elliptical, triangular, rectangular, etc. normal mode distortions 
of the jet, respectively. Propagation and growth or damping of the Fourier components can be 
described by a dispersion relation of the form 



(2) 



(1) 



Derivation of this dispersion relation is given in Appendix B. In the dispersion relation and Hn 
are Bessel and Hankel functions, the primes denote derivatives of the Bessel and Hankel functions 
with respect to their arguments. In equation (2) 
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(3b) 
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In equations (3a &; 3b) and equations (4a 8z 4b) tu|g = (w — kuj^ef and = {k — uuj^ejf?^ ■, 
7j^g = (1 — u^g/c^)~^/^ is the flow Lorentz factor, 7Aj,e = (1 ~ el^^^^^"^ Alfven Lorentz 

factor, W = p+ [r/ (r — 1)] Pjc? is the enthalpy, a is the sound speed, va is the Alfven wave speed, 
and Vms is a magnetosonic speed. The sound speed is defined by 

1/2 

VP 



a = 



f+{T&j)p/c 

where 4/3<r<5/3is the adiabatic index. The Alfven wave speed is defined by 



VA = 



Li + v;2/^2_ 



1/2 



where = BqI{A'kWq). A magnetosonic speed corresponding to the fast magnetosonic speed for 
propagation perpendicular to the magnetic field (e.g., Vlahakis & Konigl 2003) is defined by 
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3. Analytical Solutions to the Dispersion Relation 

In this section analytical solutions to the dispersion relation in the low frequency hmit, in the 
fluid and magnetic limits at resonance (maximum growth), and in the high frequency limit are 
summarized. The analytical solutions are derived in Appendix C. 

3.1. Low Frequency Limit 

Analytically each normal mode n contains a single fundamental/surface wave {u — > 0, k — > 0, 
Lo/k > 0) solution and multiple body wave {to — > 0, k > 0, to/k — > 0) solutions that satisfy the 
dispersion relation. In the low frequency limit the fundamental pinch mode (n = 0) solution is 
given by 



_ _ Uj ± Vyj 

k 1 lb VyjUj/C^ 

where the pinch fundamental mode wave speed 



msj 



+ 6 



msj 



and 



TO. 



ln( 



(5) 



(6) 



(7) 



lj J msj 

with |5| oc << 1. In this limit S is complex and this mode consists of a growing and damped 

wave pair. The imaginary part of the solution is vanishingly small in the low frequency limit. The 
above form indicates that growth, which arises from the complex value of S, will be reduced as 
(^Aj /^msj) — ^ 1- '^^^ unstable growing solution is associated with the backwards moving (in the 
jet fluid reference frame) wave. 

In the low frequency limit the surface helical, elliptical, and higher order normal modes (n > 0) 
have a solution given by 



k 



[r]Uj + Ue] ± irj 
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where rj = JjWj /^^Wg and a "surface" Alfven speed is defined by 

5? + 
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(8) 
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In equation (9) note that the Alfven Lorentz factor ^ 
n > surface wave mode perturbations when 



2 - 1 + Vl^ Jc^. Thus, the jet is stable to 



2 2 / 
Ijle [Uj 



(10) 



For example, with Uj 



c » Ue, 7e 



1> 7a,- » 7Ae 



1, B] » Bl, and using 



1 + BhATrWjCp the jet is stable when 



1 < 



1 + 
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or with Bg = Bj, We = Wj, so that vaj = VA,e, and with 7^ = 7^^, 



7A,j the jet is stable when 
(lib) 



Thus, the jet can remain stable to the surface wave modes even when the jet Lorentz factor exceeds 
the Alfven Lorentz factor. 



In the low frequency limit the real part of the body wave solutions is given by 



1/2 



X [(n + 2m - l/2)7r/2 + 



(12) 



where n specifies the normal mode, m = 1,2,3, ... specify the first, second, third, etc. body wave 
solutions, and 



en = 



Xj 13 j 



1/2 
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In the absence of a significant external magnetic field and a significant external flow e„ = as 
Xe = 7e7le^e [u^ — v\^] k"^ = 0. In this low frequency limit the body wave sohitions are either 
purely real or damped, exist only when k^^R has a positive real part, and with |e„| << 1 require 
that 



7f(u^a^)(u^-.i,) 



> . 



(13) 



Thus, the body modes can exist when the jet is supersonic and super- Alfvenic, i.e., — > 
and — > 0, or in a limited velocity range given approximately by Oj > Uj > [y^j/ (1 + 75j)]aj 
when v\j fa Oj, where jsj = (1 ~ Oj/c^)"^^^ is a sonic Lorentz factor. 

3.2. Resonance 

With the exception of the pinch fundamental mode which can have a relatively broad plateau 
in the growth rate, all body modes, and all surface modes can have a distinct maximum in the 
growth rate at some resonant frequency. 

The resonance condition can be evaluated analytically in either the fluid limit where a » Va 
or in the magnetic limit where Va » a. Note that in the magnetic limit, magnetic pressure 
balance implies that Bj = B^. In these cases a necesary condition for resonance is that 



> 



"'J 



1 - UjUe/c^ 1 + VyjjVyje/c^ 



(14) 



where Vyjj = {aj,VAj) and v^^ = (o-e^VAe) in the fluid or magnetic limits, respectively. When this 
condition is satisfied it can be shown that the wave speed at resonance is 



where 7^ = (1 — f^/c^)~^/^ is the sonic or Alfvenic Lorentz factor accompanying 
and Vwe = {cie,VAe) in the fluid or magnetic limits, respectively. 



(15) 
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The resonant wave speed and maximum growth rate occur at a frequency given by 

(2n + l)7r/4 + m7r 



(1 - Ue/V^,) - {Vyje/vl, - UeV^e/c^) 



.2^2 



1/2 ■ 



(16) 



In equation (16) n specifies the normal mode, m = specifies the surface wave, and m > 1 specifies 
the body waves. In the limit of insignificant sheath flow, = 0, and using eq. (15) for v"^ in eq. 
(16) allows the resonant frequency to be written as 



(2n + l)7r/4 + mTr 



1/2 ' 



and this predicts a resonant frequency that is primarily a function of the sound and Alfven wave 
speeds in the sheath. The effect of sheath flow is best illustrated by assuming comparable conditions 
in the spine and sheath, jwjVwj ~ IweUwe, and assuming that 'yjUj » 'yeUg in which case 



(2n + l)7r/4 + mir 



7e 



l-2{ujuj){l-vljc^)-{vl,-ui)/u) 



1/2 



The term u^/uj in the denominator indicates that the resonant frequency increases as the shear 
speed, Uj — u^, declines. In the limit 



U-i — Up 



1 - UjUe/c^ 1 + VyjjVyje/c^ 

the resonant frequency uj^^R/vwe — > oo- 

The resonant wavelength is given by A «i A*^ = 27ru^/a;*^ and can be calculated from 
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(17) 



Equations (15 - 17) provide the proper functional dependence of the resonant wave speed, frequency 
and wavelength provided (ue/uj)^ << 1 and {v^^/uj)'^ << 1. 

With the exception of the n = 0, m = 0, fundamental pinch mode, a maximum spatial growth 
rate, kY^-, is approximated by 



kf^R ^k}R = - 



1 V. 



WJ 



2 -fjUj 



ln|7^| , 



where 



\n\ 



4 iuj*^R/v^ef (1 - 2ue/uj) + (In \n\ /2y 
(ln|7^| /2)2 



1/2 



(18) 



(19) 



Equations (18) and (19) show that the maximum growth rate is primarily a function of the jet 
sound, Alfven and flow speed through v^^j/jjUj, and secondarily a function of the sheath sound, 

Alfven and flow speed through {u^^R/vwe)'^ (1 — 2ue/uj). 
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I can illustrate the dependencies of the maximum growth rate on sound, Alfven and flow speeds 
by using 



(1 - 2Ue/Uj) 



(1 - 2Ue/Uj) 



l-2{Ue/uj){l-vlJc^)-{vl,-ul)/u 



X [(2n + l)7r/4 + m7r]' 



and if say Ug = 0, then 



I7^l 



Ij ln|7^| 



1 



1 - vlJuj 



1/2 



X [(2n + l)7r/4 + m7r] 



(20) 



Thus, \Tl\ increases as increases for higher order modes with larger n and larger m and this 
result indicates an increase in the growth rate for larger n and larger m. When the sound or Alfven 
wave speed, Vyje, increases \Tl\ increases. This result indicates an increase in the growth rate at 
the higher resonant frequency accompanying an increase in the sound or Alfven wave speed in the 
sheath. 

The behavior of the maximum growth rate as the shear speed, Uj—Ue, declines is best illustrated 
by considering the effect of an increasing wind speed where {v^^ — Ug)/u'j « 1 is ignored. In this 



case 



1/2 



1 ln|7^| 4[(2n + l)7r/4 + m7r] 



(211 



and \7l\ will remain relatively independent of a;*„ even as Lo^^y^ — > oo as the shear speed decreases. 
This result indicates a relatively constant resonant growth rate as the shear speed decreases. 

In the fluid limit decline in the shear speed ultimately results in a decrease in the growth rate 
and increase in the spatial growth length. This decline in the growth rate is also indicated by 
equation (8) which, in the fluid limit, becomes 



OJ 

k 



rjUj + Ue . . ri 



1/2 



l + rj 



1 + 7? 



{Uj - Ue) . 



(22) 



Equation (22) applies to frequencies below the resonant frequency oj^m directly reveals the 
decline in growth rates as Uj — Ue — ^ 0. 

In the magnetic limit the resonant frequency lo^^R/vac — ^ oo as 



VAj + VAe 



1 - UjUe/c^ 1 + VAjVAe/c^ 

Here equation (8) indicates that the jet is stable when 

2 2/ ^2 ^ T^2 

and the jet will be stable as w*^ — > oo when 



(23) 
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,2\2 
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{l + VAjVAe/c ) 



(24) 



where I have used an equality in equation (23) in equation (8) to obtain equation (24). Equation 
(24) indicates that a high jet speed relative to the Alfven wave speed is necessary for instability. 
For example, if va = VAj = VAe and = 0, the jet is stable at high frequencies provided 



7|<(l + 4/c^)M 



(25a) 
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This high frequency condition is shghtly different from the low frequency stabihzation condition 
found when va = vaj = vac and = from equation (lib) 

7|(n,/c)2 < 471(71 - 1) . (25b) 

Note that eqs. (25a & 25b) are identical in the large Lorentz factor limit. Equations (25) predict 
that stabilization at high frequencies occurs at somewhat higher jet speeds than stabilization at 
lower frequencies. Determination of stabilization at intermediate frequencies requires numerical 
solution of the dispersion relation. A non-negligable postive external flow requires even higher jet 
speeds for the jet to be unstable. Thus, a strongly magnetized relativistic trans-Alfvenic jet is 
predicted to be KH stable and a super-Alfvenic jet can be KH stable. 



3.3. High Frequency Limit 

Provided the condition, eq. (14), for resonance is met, the real part of the solutions to the 
dispersion relation in the high frequency limit for fundamental, surface, and body modes is given 

by 

k l±VwjUj/c^ 

and describes sound waves Vy^j = aj or Alfvcn waves Vyjj = vaj propagating with and against the 
jet flow inside the jet. Unstable growing solutions are associated with the backwards moving (in 
the jet fluid reference frame) wave but the growth rate is vanishingly small in this limit. 



4. Numerical Solution of the Dispersion Relation 

The detailed behavior of solutions within an order of magnitude of the resonant frequency 
and for comparable sound and Alfven wave speeds must be investigated by numerical solution of 
the dispersion relation. Analytical solutions found in the previous section can be used for initial 
estimates and to provide the functional behavior of solutions. Numerical solution of the dispersion 
relation also allows a determination of the accuracy and applicability of the analytical expressions 
in §3. 

In this section pinch fundamental, helical surface and elliptical surface, and the associated 
first body modes are investigated in the fluid, magnetic and magnetosonic regimes. These modes 
are chosen as they have been identified with structure seen in relativistic hydrodynamic (RHD) 
numerical simulations or tentatively identified with structures in resolved AGN jets. For example, 
trailing shocks in a numerical simulation (Agudo et al. 2002) and in the 3C 120 jet (Gomez et al. 
2001) have been identified with the first pinch body mode. The development of large scale helical 
twisting of jets has been attributed to or may be associated with growth of the helical surface mode, 
e.g., 3C449 (Hardee 1981) and Cygnus A (Hardee 1996) Additionally, the development of twisted 
filamentary structures has been attributed to helical and elliptical surface and first body modes, 
e.g., 3C273 (Lobanov & Zensus 2001), M87 (Lobanov, Hardee & Eilek 2003), 3C 120 (Hardee, 
Walker &: Gomez 2005), and have been studied in RHD numerical simulations, e.g., Hardee &: 
Hughes (2003); Perucho et al. (2006). 



4.1. Fluid Limit 

In this section the basic behavior of the pinch (F) fundamental, helical (S) surface and elliptical 
(S) surface modes is investigated: (1) as a function of varying sound speed in the external sheath 
or jet spine for a fixed sound speed in the jet spine or external sheath and no sheath flow, (2) as a 
function of equal sound speeds in the jet spine and external sheath for no sheath flow, and (3) as 
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a function of sheath flow for a relatively high sound speed equal in jet spine and external sheath. 
In general only growing solutions are shown and complexities associated with multiple crossing 
solutions are not shown. For all solutions shown the jet spine Lorentz factor and speed are set 
to 7 = 2.5 and uj = 0.9165 c. Sound speeds arc input directly with the only constant being the 
sheath number density. Total pressure and spine density are quantities computed for the specified 
sound speeds. The adiabatic index is chosen to be F = 13/9 when 0.1 < aj^e/c < 0.5 consistent 
with relativistically hot electrons and cold protons (Synge 1957). For sound speeds aj^^ ~ c/VS the 
adiabatic index is set to F = 4/3. Solutions shown assume zero magnetic field. Test calculations 
with magnetic fields giving magnetic pressures a few percent of the gas pressure and Alfven wave 
speeds an order of magnitude less than the sound speeds give almost identical results. 

In Figure 1 solutions in the left column are for a fixed jet spine sound speed aj = 0.3 c and in 
the right column are for a fixed external sheath sound speed ag = 0.3 c. The solutions shown in 
Figure 1 confirm the accuracy of the low frequency solutions to the pinch fundamental mode, eqs. 
(5 & 6), and the helical and elliptical surface modes, cq. (8). Note that fast or slow wave speeds 
are possible at low frequencies depending on whether rj ~ {■jjae/jeO'j)'^ in eq. (8) is much greater 
or much less than one, respectively. The numerical solutions to the dispersion relation show that 
the maximum growth rate is primarily a function of the jet spine sound speed and only secondarily 
a function of the external sheath sound speed as indicated by eqs. (18 - 20). Where a distinct 
supersonic resonance exists, the resonant frequency is primarily a function of the external sheath 
sound speed as predicted from eq. (16). The analytical expression for the resonant frequency for the 
helical and elliptical surface modes provides the correct functional variation to within a constant 
multiplier provided Og < c/a/3 and aj < c/S. A dramatic increase in the resonant frequency and 
modest increase in the growth rate for larger jet spine sound speeds indicates the transition to 
transonic behavior. Equation (15) for the resonant wave speed and equation (17) for the resonant 
wavelength also provide a reasonable approximation to the functional variations provided < c/\/3 
and aj < c/3. These results confirm the resonant solutions found in §3.2. At frequencies more than 
an order of magnitude above resonance the growth rate is greatly reduced and solutions approach 
the high frequency limiting form given by eq. (26). Note that eq. (26) allows only relatively high 
wave speeds at high frequencies because aj < c/\/3. 

In Figure 2 the behavior of solutions to the fundamental/surface (left column) and associated 
first body mode (right column) shows how solutions change as the sound speed increases in both the 
jet spine and external sheath. Here I illustrate the transition from supersonic to transonic behavior 
for no flow in the sheath. At low frequencies the modes behave as predicted by the analytic 
solutions given in §3.1. The solutions show the expected shift to a higher resonant frequency that 
is primarily a function of the increased external sheath sound speed and an accompanying increase 
in the resonant growth rate that is primarily a function of the increased jet spine sound speed. The 
resonance disappears as sound speeds approach c/^/3 as the jet becomes transonic as predicted 
by the resonance condition in §3.2. In the transonic regime high frequency fundamental/surface 
mode growth rates and wave speeds are identical with wave speeds given by eq. (26). Provided the 
jet is sufficiently supersonic, i.e., aj^e < 0.5 c, the maximum growth rate of the first body mode 
is greater than that of the pinch fundamental mode, is comparable to that of the helical surface 
mode, and is less than that of the elliptical surface mode. A narrow damping peak shown for the 
helical first body (Bl) solution when aj^e = 0.4 c is indicative of complexities in the body mode 
solution structure. In the transonic regime growth of the first body mode is less than that of the 
pinch fundamental, helical surface and elliptical surface modes. 

Figure 3 illustrates the behavior of fundamental/surface and first body modes as a function 
of the sheath speed for equal sound speeds in spine and sheath of aj^e = 0.4 c. For this value 
of the sound speeds a sheath speed = provides a supersonic solution structure baseline. At 
low frequencies the surface modes behave as predicted by eq. (8), and the wave speed rises as Ug 
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Fig. 1. — Solutions to the dispersion relation for pinch fiindamental, helical surface and elliptical surface modes 
for different sound speeds in the sheath (left column) and in the spine (right column) are shown for no sheath flow. 
The real part of the wavenumber, krRj, is shown by the dashed lines and the imaginary part , kiRj, is shown by the 
dash-dot lines as a function of the dimensionless angular frequency, uiRj/uj. For the pinch mode the vertical lines 
indicate the maximum growth rate range. Otherwise, the vertical lines indicate the location of maximum growth. 
Immediately under the dispersion relation solution panel is a panel that shows the relativistic wave speed, ywVw/c. 
Line colors indicate the sound speed in units of c: (black) 0.10, (blue) 0.20, (cyan) 0.30, (green) 0.40, & (red) 0.577. 
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Fig. 2. — Solutions to the dispersion relation for pinch fundamental, helical surface, elliptical surface (left column) 
and the first body (right column) modes are shown for equal sound speeds in spine and sheath and no sheath flow. 
Real and imaginary parts of the wavenumber as a function of angular frequency are shown as in Figure 1. Locations 
of the maximum growth rate are indicated by the vertical solid lines. A vertical arrow (helical Bl) indicates a narrow 
damping feature. The underlying panel shows the relativistic wave speed, 'ywVw/c. Line colors indicate the sound 
speed in units of c: (black) 0.10, (blue) 0.20, (green) 0.40, & (red) 0.577. 
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Fig. 3. — Solutions to the dispersion relation for pinch fundamental, helical surface, elliptical surface (left column) 
and first body (right column) modes as a function of the sheath speed for equal sound speeds in spine and sheath. 
Real and imaginary parts of the wavenumber as a function of angular frequency are shown as in Figure 1. Locations 
of the maximum growth rate are indicated by the vertical solid lines. Vertical arrows (helical Bl) indicate damping 
features. The underlying panel shows the relativistic wave speed, jwVw/c. Line colors indicate the sheath speed in 
units of c: (black) 0.0, (blue) 0.20, (cyan) 0.35, (green) 0.40, & (red) 0.60. 
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increases. As Ue increases the resonant frequency increases in accordance with eq. (16). On the 
other hand, the growth rate at resonance does not vary significantly in accordance with eqs. (18 &: 
19). When the sheath speed exceeds the sound speed, solutions make a transition from supersonic 
to transonic structure. Note that the transition point between supersonic and transonic behavior 
is similar but not identical for the helical and elliptical surface modes, i.e., ocurs at a slightly 
lower sheath speed for the elliptical mode. The first body modes also show an increase in resonant 
frequency with little change in the maximum growth rate provided the sheath speed remains below 
the sound speed. A significant damping feature in the helical first body (Bl) panel, is found. While 
a similar damping feature was not found for the pinch and elliptical first body mode, this does not 
indicate a significant difference as the root finding technique docs not find all structure associated 
with the body modes. The body mode solution structure is complex with multiple solutions not 
shown here and modest damping or growth can occur where solutions cross, e.g., Mizuno, Hardee 
& Nishikawa (2006). When the sheath speed exceeds the sound speed the maximum body mode 
growth rate defines significantly. This result is quite different from the transonic solution behavior 
illustrated in Figure 2 when aj^e = 0.577 c for no sheath flow. Thus, sheath flow effects stability 
of the relativistic jet beyond that accompanying an increase to the maximum sound speed in the 
absence of sheath flow. The reduction in growth of the body modes in the presence of sheath flow 
provides the relativistic jet equivalent of non-relativistic transonic/subsonic jet solution behavior. 
At the higher frequencies wave speeds are identical, with wave speeds given by eq. (26). Note that 
the high frequency wave speeds are nearly independent of Ug- 

4.2. Magnetic Limit 

In this subsection the basic behavior of pinch, helical and elliptical modes is investigated: (1) 
as a function of varying Alfven speed in the external sheath or jet spine for a fixed Alfven speed 
in the jet spine or external sheath and no sheath fiow, (2) as a function of equal Alfven speeds 
in the jet spine and external sheath for no sheath fiow, and (3) as a function of sheath speed for 
a relatively high Alfven speed equal in jet spine and external sheath. In general only growing 
solutions are shown and complexities associated with multiple crossing solutions are not shown. 
For all solutions shown the jet spine Lorentz factor and speed are set to 7 = 2.5 and uj = 0.9165 c. 
Alfven speeds arc on the order of two magnitudes larger than the sound speed and arc determined 
by varying the sound speeds but with a gas pressure fraction on the order of 0.01% of the total 
pressure. Only the sheath number density is held constant. The adiabatic index is set to F = 5/3 
when aj^e/c << 0.1 consistent with low gas pressures and temperatures. 

The solutions shown in Figure 4 confirm the theoretical predictions in the magnetic limit with 
behavior depending on the Alfven speed like the behavior found for the sound speed (see Figure 
1). The pinch fundamental mode (not shown) has a growth rate almost entirely dependent on 
sound speeds and is negligable in the magnetic limit as predicted by eq. (6). In Figure 4 solutions 
in the left column are for a fixed jet spine Alfven speed VAj = 0.3 c and in the right column are 
for a fixed external sheath Alfven speed vab = 0.3 c. The solutions shown confirm the accuracy 
of the low frequency solutions for helical and elliptical surface modes given by eq. (8). Note that 
low frequency wave speeds can be high or low depending on the values of rj = 'yjWj/j^We, VAelle 
and VAjllj- The numerical solutions to the dispersion relation show that the maximum growth 
rate is primarily a function of the jet spine Alfven speed and only secondarily a function of the 
external sheath Alfven speed as predicted by eqs. (18 - 20). The resonant frequency is primarily a 
function of the external sheath Alfven speed as predicted by eq. (16). The analytical expression for 
the resonant frequency of the helical and elliptical surface modes provides the correct functional 
variation to within a constant multiplier provided VAj,e < 0.5 c. Decrease in the growth rate 
for jet sheath Alfven speeds VAe > 0.5 c indicates the transition towards trans-Alfvenic behavior. 
Equation (15) for the resonant wave speed and equation (17) for the resonant wavelength also 



- 14 - 



Helical S (v = 0.3) 




Helical S (v = 0.3) 



i^-V " ^ 


0.3 
i0,5 


0.80 

\ 

\ : 

\ , 

\ — 1 1 1 ' r' ^ 


1 1 1 1 1 \ 








1 



0.1 1 10' 

Elliptical S [y J c = 0.3) 



% - 0-10-c 


Ml 

.50 


'^^^^^ - 
0.80 
" \ 
. "\ 
. '\ 

\ 

\ : 

— \ — 1 1 1 1 ' 1 ' — ^^-^ — ■ 


1 1 1 1 1 — \ — II 

— ....I 


— 

1,1 





cjR/u. 

Fig. 4. — Solutions to the dispersion relation for pinch fundamental, helical surface and elliptical surface modes 
for different Alfven speeds in the sheath (left column) and in the spine (right column) are shown for no sheath flow. 
Sound speeds are aj,e '~ O.OluAj.e. As in Figures 1-3, the real part of the wavenumber, krRj, is shown by the 
dashed lines and the imaginary part , kiRj, is shown by the dash-dot lines as a function of the dimensionless angular 
frequency, uiRj/uj. The vertical lines indicate the location of maximum growth. The underlying panel shows the 
relativistic wave speed, "/mV-w/c. Line colors indicate the Alfven speed in units of c: (black) 0.10, (cyan) 0.30, (green) 
0.50, & (red) 0.80. 



provide a reasonable approximation to the functional variations for v^j^e < 0.5 c. At frequencies 
more than an order of magnitude above resonance the growth rate is greatly reduced and solutions 
approach the high frequency limiting form given by eq. (26). The surface modes have relatively 
slow wave speeds, 'jwVw/c < 1 at high frequencies when the Alfven wave speed vaj > 0.5 c. Unlike 
the fluid case, the helical and elliptical surface modes are stabilized for Alfven speeds somewhat in 
excess of VAj,e ~ 0.8 c in accordance with eqs. (8 & 24). 

In Figure 5 the behavior of solutions to the pinch fundamental mode is shown in addition 
to the helical and elliptical surface (left column) and associated first body modes (right column) 
and the figure shows how solutions change as the Alfven speed increases in both the jet spine and 
external sheath. The sound speed is aj^e = 0.2 c for the pinch fundamental mode panel in order to 
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Fig. 5. — Solutions to the dispersion relation for pinch fundamental, helical surface, elliptical surface (left column) 
and the first body (right column) modes are shown for equal sound speeds in jet and sheath and no sheath flow. 
Pinch fundamental mode sound speed is a^.e = 0.2 c. Sound speeds for all other cases are a^^e ~ O.Olv Aj,e- Real and 
imaginary parts of the wavenumber as a function of angular frequency are shown as in Figure 4. Locations of the 
maximum growth rate are indicated by the vertical solid lines. A vertical arrow (elliptical Bl) indicates a narrow 
damping feature. The underlying panel shows the relativistic wave speed, 'y-wVw/c. Line colors indicate the Alfven 
speed in units of c: (black) 0.10, (blue) 0.20, (green) 0.40, & (red) 0.60. 
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Fig. 6. — Solutions to the dispersion relation for pinch fundamental, helical surface, elliptical surface (left column) 
and the first body (right column) modes are shown for equal sound speeds in jet and sheath for different sheath 
flow speeds. The pinch fundamental mode sound speed is aj,e — 0.2 c. Sound speeds for all other cases are 
a.j,e ~ O.OlwAj.e ~ 0.005 c. Real and imaginary parts of the wavenumber as a function of angular frequency are 
shown as in Figure 4. Locations of the maximum growth rate are indicated by the vertical solid lines. A vertical 
arrow indicates low frequency damping of the pinch Bl solutions. The underlying panel shows the relativistic wave 
speed, ywVw/c. Line colors indicate the sheath speed in units of c: (black) 0.0, (blue) 0.20, (cyan) 0.30, (green) 0.40, 
& (red) 0.60. 
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illustrate the mode behavior with increasing Alfven speed. Sound speeds for all body modes and 
for helical and elliptical surface modes are aj^e ~ O.Olf Aj,e- Here the transition from super- Alfvenic 
towards trans-Alfvenic behavior for no flow in the sheath is illustrated. At low frequencies the 
modes behave as predicted by the analytic solutions given in §3.1. The growth rate of the pinch 
fundamental mode is reduced as the Alfven speed increases as predicted by eq. (6). The surface 
and body mode solutions show the expected shift to a higher resonant frequency that is primarily 
a function of the increased sheath Alfven speed and an accompanying increase in the resonant 
growth rate that is primarily a function of the increased spine Alfven speed. The resonance moves 
to higher frequency but the maximum growth rate is reduced for Alfven speeds VAj,e > 0-60 c 
and all modes become stable at higher Alfven speeds in accordance with eqs. (8 &: 24). At high 
frequencies wave speeds are given by eq. (26). Provided the jet is sufficiently super-Alfvenic, i.e., 
'fAj,e < 0-6 c, the maximum growth rate of the first body mode is much greater than that of the 
pinch fundamental mode, is comparable to that of the helical surface mode, and is less than that 
of the elliptical surface mode. A narrow damping peak shown for the elliptical body mode (Bl) 
solution when VAj,e = 0.6 c indicated by the arrow is indicative of complexities in the body mode 
solution structure. 

Figure 6 illustrates the behavior of fundamental/surface and first body modes as a function 
of the sheath speed for an equal Alfven speed in spine and sheath of VAj,e = 0.5 c. For this value 
of the Alfven speeds a sheath speed Ug = provides a super-Alfvenic solution structure baseline. 
The sound speed is aj^e = 0.2 c for the pinch fundamental mode panel in order to illustrate the 
mode behavior with increasing sheath speed. Sound speeds for all other cases are aj^e ~ O.OlfAj.e- 
Solutions in the first pinch body mode panel show the damping solution as opposed to the purely 
real solution at the lower frequencies (indicated by the arrow) . At higher frequencies the body mode 
is growing. At low frequencies the surface modes behave as predicted by the analytic solutions given 
in §3.1 and the growth rate of the surface modes decreases as Ug increases. Additionally, the growth 
rate at resonance decreases as expected for this relatively high Alfven speed as the sheath speed 
increases. At the higher frequencies wave speeds are identical, with wave speeds given by eq. (26). 
Note that the high frequency wave speeds are relatively independent of u^- When the velocity shear 
speed drops to less than the "surface" Alfven speed, see eq. (lib), the helical and elliptical surface 
modes and the first body modes are stabilized. This surface and body mode mode stabilization 
occurs when sheath speeds exceed ~ 0.5 c. However, note that the maximum pinch fundmental 
mode growth rate is insensitive to the sheath speed and remains unstable at Ug = 0.6 c even when 
all other modes are stabilized. 

4.3. A High Sound and Alfven Speed Magnetosonic Case 

In this subsection the basic behavior of the pinch fundamental, helical surface, elliptical surface 
and associated first body modes is illustrated for different sheath speeds. The sheath speeds span a 
solution structure from supersonic to transonic but still super-Alfvenic flow. Here the sound speed 
in jet spine and external sheath are set equal with aj^c = 0.577 c and Alfven speeds are set equal 
with VAj,e = 0.5 c. The solutions for this case are shown in Figure 7. With no sheath flow the 
fundamental/surface and flrst body modes show a typical supersonic and super-Alfvenic structure 
albeit the pinch fundamental mode now has a maximum growth rate comparable to the helical 
and elliptical surface modes as a consequence of the high sound speed. The associated first body 
modes also have maximum growth rates comparable to the fundamental/surface modes. Increase 
in the sheath speed results in a decrease in the growth rate of the helical and elliptical surface 
modes at low frequencies as predicted by eq. (8). The low frequency growth rate of the pinch 
fundamental also declines with increasing sheath speed. The resonant frequency increases with 
increasing sheath speed as expected from the analytical and numerical studies performed in the 
fiuid and magnetic limits and the fundamental/surface modes take on a transonic structure for 
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Fig. 7. — Solutions to the dispersion relation for pinch fundamental, helical surface, elliptical surface (left column) , 
and the associated first body (right column) modes are shown for a maximal spine and sheath sound speed, aj^e = 
0.577 c, and a slightly smaller spine and sheath Alfven speed, VAi,e = 0.5 c, for different sheath flow speeds. As in 
previous flgures the real part of the wavenumber, krRj, is shown by the dashed lines, the imaginary part , kiRj, 
is shown by the dash-dot lines, and the vertical lines indicate the location of maximum growth. Arrows indicate 
damping features. The underlying panel shows the relativistic wave speed, 'ywVm/c. Line colors indicate the sheath 
speed in units of c: (black) 0.0, (blue) 0.20, (green) 0.40, & (red) 0.50. Fast and slow refer to the faster and slower 
moving solutions and the yellow extension indicates a damped solution. 
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sheath speeds 0.4 c < Ue < 0.1 c. 

At high frequencies the fundamental/surface modes exhibit very high growth rates provided 
sheath flow remains below the Alfven speed. On the other hand, the maximum growth rate of the 
first body modes declines as the sheath speed increases and is reduced severely when Ue > 0.1 c. 
This behavior is similar to what is found for non-relativistic jets as flow enters the transonic and 
super- Alfvenic regime (Hardee &: Rosen 1999). Additional increase in the sheath flow speed to 
Ue > 0.4 c results in a decrease in the growth rate of the fundamental/surface modes. Solutions 
for the helical and elliptical surface modes shown in Figure 7 for a sheath speed = 0.5 c equal 
to the Alfven speed illustrate some of the complexity associated with barely super-Alfvenic flow. 
Here limited growth is associated with both the slow and fast helical and elliptical surface solution 
pair. At slower sheath speeds in the super-Alfvenic regime growth is associated with the slow 
surface solution, i.e., backwards moving in the jet fluid reference frame. The yellow dash-dot 
line extension at higher frequencies in the helical and elliptical surface panels indicates a damped 
solution. Solutions were very difficult to follow in this parameter regime and it is possible that 
some solutions were not found. When the sheath speed Ue > 0.5 c all modes are stabilized. 

A choice of Alfven speeds greater than sound speeds results in a more magnetic like solution 
structure like that shown in §4.2. A choice of Alfven speeds more than a factor of two less than sound 
speeds produces a more fluid like solution structure like that shown in §4.1. The more complicated 
solution structure illustrated in Figure 7 only occurs for a relatively narrow range of high sound 
speeds with similar or slightly lesser Alfven speeds. In general, the detailed solution structure for 
situations in which sound and Alfven speeds are comparable must by examined individually, e.g., 
Mizuno, Hardee &; Nishikawa (2006), and further investigation of these cases is beyond the scope 
of the present paper. 



5. Summary 

The analytical and numerical work performed here provides for the first time a detailed analysis 
of the KH stability properities of a RMHD jet spine-sheath configuration that allows for relativistic 
motions of the sheath, sound speeds up to c/\/3, and, by keeping the displacement current in the 
analysis, Alfven wave speeds approaching lightspeed and large Alfven Lorentz factors. In the fluid 
limit, the present results confirm an earlier more restricted low frequency analytical and numerical 
simulation study performed by Hardee & Hughes (2003). Provided the jet spine is super-sonic and 
super-Alfvenic internally and also relative to the sheath, the helical, elliptical and higher order 
surface modes and the pinch, helical, elliptical and higher order first body modes have a maximum 
growth rate at a resonant frequency. The pinch fundamental growth rate is significant only when 
the sound speeds, Uj^e ~ c/\/3. In general, the first body mode maximum growth rate is: greater 
than the pinch fundamental mode, slightly greater than the helical surface mode, slightly less than 
the elliptical surface mode, and occurs at a higher frequency than the maximum growth rate for 
the fundamental/surface mode. 

The basic KH stability behavior as a function of spine-sheath parameters is indicated by the 
analytic low frequency surface mode solution and by the behavior of the resonant frequency. The 
analytic surface mode solution valid at frequencies below resonance is given by 



k k k (^i + vl/j^ec'^)+v{l + Vl/j]c'^ 

where 



1/2 

— (27) 



^l-(7i.^.+7ie^e)|^, (28) 
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and ri = T^VFj/^/e I'Fe , Vl = B'^/AtiW, W = p + [r/ (F - 1)] P/c^ and 7^ = (1 - vyc^)-^/'^. 
Equation (27) provides a temporal growth rate, uJi{k), and a wave speed, Vw = ujr/k. The reciprocal 
provides a spatial growth rate ki{u}), and growth length £ = k^^. Increase or decrease of the 
growth rate, dependence on physical parameters and stabilization at frequencies/wavenumbers 
below resonance is directly revealed by LOi in eq. (27). Note that higher jet Lorentz factors reduce 
LOi through the dependence on 77. 
The resonant frequency is 

^ 0^ 7 TTT^ ' (29) 



(1 - Ue/v^ f - {Vu,e/v^ - UeVuie/c^Y 



where v^^ is the wave speed at resonance, eq. (15). The resonant frequency increases as the sheath 
sound or Alfven wave speed, v^e = {o,e,VAe) increases and lo* — ^ 00 when the denominator 
decreases to zero as 

Uj tig Vwj ~t~ '^we 

where Vyjj^e = (oj,e) VAj,e) in the fluid and magnetic limits, respectively. Since eq. (27) applies below 
resonance the overall behavior of the growth rate is indicated by Ui. Thus, growth rates decline to 
zero as {uj — lig)^ — ^s/7j7e — ^ 0- numerical analysis of the dispersion relation shows that 
the pinch fundamental and all first body modes are comparably or more readily stabilized and thus 
the jet is KH stable when 

{uj-u,f-Vl,h]^l<Q. (30) 

This stability condition takes on a particularly simple form when conditions in spine and sheath 
are equal, i.e., = Bj, We = Wj, so that vaj = VA,e, and with 7^ = 7A,e = 7Aj 

7|7e - Uef < 47i ill - 1) (31) 

indicates stability. This result implies that a trans- Alfvenic relativistic jet with jjuj > Java will 
be KH stable, and that even a super- Alfvenic jet with jj » ja can be KH stable. 



6. Discussion 

Formally, the present results and expressions apply only to magnetic fields parallel to an axial 
spine-sheath flow in which conditions within the spine and within the sheath are independent of 
radius and the sheath extends to infinity. A rapid decline in perturbation amplitudes in the sheath 
as a function of radius, governed by the Hankel function in the dispersion relation, suggests that 
the present results will apply to sheaths more than about three times the spine radius in thickness. 

The relativistic jet is transonic in the absence of sheath flow only for spine and sheath sound 
speeds ~ c/\/3. Only in this regime does the pinch fundamental have a significant growth rate and, 
in general, we do not expect the pinch fundamental to grow significantly on relativistic jets. On 
the other hand, the pinch first body mode can have a significant maximum growth rate and would 
dominate any axisymmctric structure. The elliptical and higher order surface modes have increas- 
ingly larger maximum growth rates at resonant frequencies higher than the helical surface mode, 
and the maximum first body mode growth rates for helical and elliptical modes are comparable to 
that of the surface modes. Nevertheless, we expect the helical surface mode to achieve the largest 
amplitudes in the non-linear limit as a result of the reduced saturation amplitudes that accompany 
the higher resonant frequency and shorter resonant wavelengths associated with the higher order 
surface modes and all body modes. 
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In astrophysical jets we expect a toroidal magnetic field component, and possibly an ordered 
helical structure and accompanying flow helicity. Jet rotation (e.g., Bodo et al. 1996), or a radial 
velocity profile (e.g., Birkinshaw 1991) will modify the present results but will not stabilize the 
helical mode. Two dimensional non-rclativistic slab jet theoretical results, indicate that KH sta- 
bilization occurs when the velocity shear projected on the wavevector is less than the projected 
Alfven speed (Hardee et al. 1992). In the work presented here magnetic and flow field are parallel 
and project equally on the wavevector which for the helical (n = 1) and elliptical (n = 2) mode 
lies at an angle 9 = tan~^{n/kR) relative to the jet axis. Provided magnetic and flow helicity and 
radial gradients in jet spine/sheath properties are not too large we expect the present results to 
remain valid where itj g and Bj g refer to the poloidal velocity and field components. 

KH driven normal mode structures move at less than the jet speed. The fundamental pinch 
mode moves backwards in the jet frame at about the sound speed nearly independent of the sheath 
properites and thus moves at nearly the jet speed in the source/observer frame. Low frequency and 
long wavelength helical and higer order surface modes are advectcd with wave speed indicated by eq. 
(27) and move slowly in the source/observer frame for light, i.e. 7] = ^"jWj /^IW^ < 1, and/or for 
magnetically dominated flows. Higher frequency (above resonance) and shorter wavelength normal 
mode structures move backwards in the jet frame at the sound/ Alfven wave speed, have a wave 
speed nearly independent of the sheath properties, and can move slowly in the source/observer 
frame only for magnetically dominated flows. 

Where flow and magentic fields arc parallel, current driven (CD) modes are stable (Isotomin 
&L Pariev 1994, 1996). Where magnetic and flow fields are helical CD modes can be unstable 
(Lyubarskii 1999) in addition to the KH modes. CD and KH instability are expected to produce 
helically twisted structure. However, the conditions for instability, the radial structure, the growth 
rate and the pattern motions are different. For example, KH modes grow more rapidly when the 
magnetic field is force-free (e.g., Appl 1996), and non-relativistic simulation work (e.g., Lery et 
al. 2000; Baty & Keppens 2003; Nakamura & Meier 2004) indicates that CD driven structure is 
internal to any spine-sheath interface and moves at the jet speed. 

The differences between KH and CD instability can serve to identify the source of helical 
structure on rclativistic jets and allow determination of jet properties near to the central engine. 
Perhaps the observation of relatively low proper motions in the TeV BL Lacs when intensity mod- 
eling requires high flow Lorentz factors (Ghisellini et al. 2005) is an indication of a magnetically 
dominated KH unstable spine-sheath configuration. 
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A. Linearization of the RMHD Equations 

In vector notation the relativistic MHD continuity equation, energy equation, and momentum 
equation can be written as: 

[7/9]+V-[7Pv]=0, (Al) 
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dt 
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These equations along with Maxwell's equations 



V • B = V • E = Airpq 

VxB = i^E + ^j VxE = -i^] 
c at c •' c at 

and assuming ideal MHD with comoving electric field equal to zero 



V X B 
E = , 

c 

provide the complete set of ideal RMHD equations. In the above W is the enthalpy, the Lorentz 
factor 7 = (1 — V • v/c^)~^/^, and p is the proper density. In what follows I will assume that the 
efi^ects of radiation can be ignored, the enthalpy is given by 

r P 

and the condition for isentropic flow is given by 



5^ + --^ 



The general approach to analyzing the time dependent properties of this system is to linearize 
the ideal RMHD equations, where the density, velocity, pressure and magnetic field are written as 
P = + Pi ; V = u + vi (we use vq = u for notational reasons) , P = Pq + Pi E = Eq + Ei , and 
B = Bo + Bi, where subscript 1 refers to a perturbation to the equilibrium quantity with subscript 
0. Additionally, W = Wq + Wi, j'^ = (70 + 71)^ — 7o + 27oU • vi/c^ and 71 ~ 70 u • Vi/c^. It 
is assumed that the initial equilibrium system satisfies the zero order equations. The linearized 
continuity, energy and momentum equation become 

d 

^ [loPi + IlPo] + V • [70P1U + 70P0V1 + 7iPou] = , (A4) 

I [j^Wi - Pi/c^ + 27^ (u • vi/c2) Wo] + V • [ilWiVi + 274 (u • vi/c2) Won + 7o'W^oVi] 

[^l • vi/c') + (1 + ^Vc')Bo-Bi - (u • Bi/c + vrBo/c)u • Bo/c] (A5) 
■ [2(Bo-Bi)u+B2vi - (u • Bo)Bi - (u • Bi)Bo - (vrBo)Bo] = , 



and 



The linearized Maxwell equations become: 

V • Bi = V • El = Airpqi 

VxBi = l|Ei + fji VxEi = -i|i 
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where I keep the displacement current in order to allow for strong magnetic fields and Alfvcn wave 
speeds comparable to lightspeed. Under the assumption of ideal MHD, the comoving electric field 
is zero, the equilibrium charge density pqfi = 0, and the electric field 



El 



u X Bi + vixBo 



is first order, the charge density pqi = (v • Ei) /47r is also first order, and the electrostatic force 
term, pgiEi, is second order and dropped from the linearized momentum equation. The condition 
for isentropic perturbations becomes 



Pi 



Po 



Pi 



This basic set of linearized RMHD equations is similar to those found in Bcgclman (1998) but 
allows a relativistic zeroth order velocity, i.e., v = u + vi and u < c whereas Begelman allowed 
only for relativistic first order motions, Vi. 

In what follows let us model a jet as a cylinder of radius R, having a uniform proper density, 
Pj, a uniform axial magnetic field, Bj = -B^j-, and a uniform velocity, uj = u^j. The external 
medium is assumed to have a uniform proper density, pg, a uniform axial magnetic field, Bg = B^^ei 
and a uniform velocity, Ug = u^^e- An external velocity could be the result of a wind or sheath 
outflow around a central jet, Ug > 0, or could represent backflow, Ue < 0, in a cocoon surrounding 
the jet. The jet is established in static total pressure balance with the external medium where the 
total static uniform pressure is P* = Pe + B^/Stt = Pj = Pj + Bj/Sir. Under these assumptions 
the linearized continuity equation becomes 



^ boPi + liPo] + u-^ [70/91 + 7iPo] + ToPoV • vi = . 



The linearized energy equation becomes 



d_ 

dt 



70^1 



P 



1 , o 4.UVzl„. 

^2+ 270^^0 



+ u 



d_ 

dz 



+ -fiWoV • vi = . 



(A7) 



(A8) 



This result for the linearized energy equation is found by noting that the magnetic terms in the 
energy equation linearize to 

Bo + uf^B.i] - {uBo)V • Bi + . vi - 

-^0 [V • VI - + 52 [V • VI - £v,,] 

where I have used 



= 



—5^1 + u—B^i = 

ot oz r 



d . . d ' 



-Bn 



V • vi - — u^i 
oz 



from d^i/dt = Vx(uxBi) + Vx (vixBq). The linearized momentum equation becomes 



at 



vi + u • Vvi — 



1 



(|vixBo) xBo 



-VA - ^l^'i + ^ [(V X Bo) X Bi + (V X Bi) X Bo] + 3^ (ux|Bi) xBo 



(A9) 
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where I have used 



and 



jpxBi ^ (V X Bo)xBi 
c 47r 



jixBo (VxBi)xBo 1 „ (VxBi)xBo^ 1 
— -Jill xt$o = : h 



Ancdt 



47r 



u X — Bi + — vixBo ) xBo , 



which includes the displacement current. The components of the linearized momentum equation 
can be written as 



1 + 



Vi 



76 



d 



/ dt 



d 
oz 



rl 



d_ 
dr 



Pi + 



Bo f d 



4tt V dz 



-Brl + 



u d 



Brl 



d_ 
dr 



B 



zl 



(AlOa) 



76 



d 



d 



1 d 



Pi + 



Bo f d 



Att \ dz 



B,pi + 



u d 



B,. 



1 d 



and 



d a 



dz (? dt 



(AlOb) 
(AlOc) 



where = BI/{A-kWo). 



B. Normal Mode Dispersion Relation 

In cylindrical geometry perturbations pi, vi, Pi, and Bi can be considered to consist of Fourier 
components of the form 

/i(r,0,z,t)=/i(r)e*(^-±«^-*) 

where the flow is in the z direction, and r is in the radial direction with the jet bounded by r = i?. 
In cylindrical geometry n, an integer, is the azimuthal wavenumber, for n > waves propagate at 
an angle to the flow direction, where +n and —n refer to wave propagation in the clockwise and 
counterclockwise sense, respectively, when viewed outwards along the flow direction. In general the 
goal is to write a differential equation for the radial dependence of the total pressure perturbation 
PI = Pi + (Bi • Bo)/47r = P*{r) exp[i{kz zt ncj) — tot)]. The differential equation can be obtained 
from the energy equation by using the momentum equation and writing the velocity components 
VriiV(j,i-,Vzi in terms of P*,u,Bq. The components of the linearized momentum equation (eqs. 
AlOa, b, k. c) written in the form 



and 



70^0 



1 + 



1 + 



d 



d 



d 



d 



d 



d 



dt dz 



dr 
1 d 



Att \ dz 



-P* jf- 
r dcj) ^ Ait 



u d 



dt 



dt 



Bo f d 



-Ba)1 + 



d 



dz"t"^7^dt^t" 



Bo 



d 



u d 



Att \ dz 



dt 



(Bla) 
(Bib) 

(Blc) 



along with 





— Bi = -c(V X El) = Vx(u X Bi) + Vx(vixBo) 
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are used to provide relations between Bri and Vri, and B^i and v^pi 

d d d 

— iJ^l + U—Brl — DQ—Vrl , 

and to provide a relation between B^i, Vzi, and P* 

iBq f 1 , r 1 \ r 9 p* I 9 p*] _ -Bp 9 p* 
To obtain equation (B2c) I have used 



— + u—B^i = 

ot az r 



where 



-7o^oV • vi = 7o + 



1 



r 1 



^a^ r-lc2 
from the energy equation (eq. A8), and 

Wi=pi + 



d , . d ' 



^ dz ^ dt 



-Bn 



d 



Pi + 270^4^0 



(B2a) 
(B2b) 



(B2c) 



d d 



r Pi 



r - 1 c2 

Using equations (Bla, b, 8z c) combined with 



1 r 1 



§lfi{r,(t),z,t) = -za;/i(r)e^('=^±"'^-^*) 
|:/i(r,</.,z,t) = |:/i(r)e^('=^±"'^-*) 

allows the velocity components to be written as 



ku-oj{l+ ^ 



ku- + 



^ P* ^ .^0 /, u\ 



Brl 



and 



7o'H^o - a;] v,^ = - (k - a; J) P^* + ^ - ^zi 



P, 



01 



(B3a) 
(B3b) 

(B3c) 
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The perturbed magnetic field components from equations (B2a, b,&; c) become 

kVrl 



B, 



rl 



ku — LO 



-Bn 



B. 



kv, 



01 



ku — LO 



Bo , 



and 



B 



zl 



1 

Wo 


'i^-{k-uu/c^) ^ 


Pt + ll[{k- 


uu/c?) + {ku — u) ^] Vzi 


{ku -io) + Vj 


Uku^ 







(B4a) 
(B4b) 

So (B4c) 



wfiere I fiave used 



k + {ku — u) u/c^ = 7o [{k — uu/cP) + {ku — u) u/cP] , 

7^ {ku - uj) + r (r - 1)"^ c-2j + oj/c^ = 7^ [{ku - uj) ja^ - {k - uju/c^) u/c^ 



and 



a = ^ + 



r 1 1 



VP 



to obtain the expression for B^i- Using equations (B4a, b, & c) for the perturbed magnetic field 
components, I obtain the following relations between the perturbed velocity components Vi and 
the total pressure perturbation P^: 



{ku — a;) 



Vri = Cr-i-Pi = i^-i-Pi = i r- 

dr Xdr j^WojI {ku - uf - {k - uu / c^f v\ 



{ku — u) 



9 p* 



v^^^C^P^ = T-^P^ = T- r 

''llWo-il {ku-uf -{k-uu/c^fv\ 



and 



= C,Pl = - 



{ku — oj){k — (jJu/(?) 



llWo{{ku - uf + 7>i -{k- uu/c^f] } 

To obtain the above relationships I have used 



1 ) 



p* 
^1 ■ 



(B5a) 
(B5b) 

(B5c) 



/ k — oju/c^ oj 



{ku — oj) 4-1 

7q \ ku — LO 



+ = 7I 



{k-uju/c^f 2 

(ku — uj) r Va 

^ ' {ku-uj) ^ 



in addition to 

7^ {ku - w) \ar'^ + r (r - l)"^ c-^j + w/c^ = 7^ {{ku - u)) jo? - {k- tou/c^) u/c^] 



where 



1 + Vj/c2 
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is the Alfven wave speed and 7^ = (l — ''^a/c^) ^ is an Alfvenic Lorentz factor. Note that V\ = 
^\v\ and 7^ = 1 + V^/c^. Thus we have that 



r* 9^ -p* \ Cr 9 p* _i_ d p* _i_ 



i 8^ p* I 1 i 9 p* i p* i p* 



(B6) 



Using the energy equation (eq. A8) written in the form 



inserting 



d d 

— + u—Bzi = —Bq 
at oz 



d 
oz 



and using Vzi = C^P^* gives 



— z 



1 + 



V • vi = 

{ku-uj) 



7^ (feu-w) 



27q (feti - 



%Y-. ' 



To (fcw-'^). 



(B7) 



where 



(B8) 



Setting equations (B6) and (B7) equal gives us a differential equation for Pj* in the form of Bessel's 
equation 

92 d 



a;:2^i+-^A* + [/3V-n2]pr = o 



dr 



(B9) 



where 



+kXCz + 



VI Y 

7(f (ku-ui) 



1 + ^. y 



To 



To 



I can simplify the expression for /J^ by writing 



(BIO) 



X ({ku- co) 

{ku-uj) + ^Y\ ^oWo 

To 



Y + 



iku-u) + ^Y 
7o . 



ka + 



270 - ^ — 



y2 

To 



from which it follows that 
X 



(ku + ^Y 

To 



(ku — uj) 

' ^ Y + 7o (A;tt - w) [(fc - wn/c^) + {ku - cv) u/c^] 
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where I have used 27q {ku — lo) ujc? = 7o [(^ — uu/c?) + {ku — u) u/c^] — k. Substituting the ex- 
pressions for X and from equations (B5a, b, & c), and Y from equation (B8) and modest 
algebraic manipulation yields 



(a'^+'y\v\){ku—u>)'^+'^\v\a'^(ku-i^)(k-Luu/c^)u/c'^ 



{ku — oj)^ + {ku — to) {k — uu/cP^ a^u/cP 



(ku—u>)^ \^{k~LL>u/ c^'j^ a?- +(ku—uj){k—uiu/ c^^ua? / c 
{a?WA^A)i^u-^f-lA'"'Wik-U}u/c'^f 



(Bll) 



Additional regrouping provides the following form 



1^ {a^ +'y\v\^{ku—LL))'' +^\v\a'^ {ku—uj){k~u)u/ c^)u/ j 
[{a'^+f\v\){ku-wf+^\v\a^{ku-uj){k-uju/c^)u/c'^] \^{ku-wf-{k-wu/c'^)'^a^'^ 



(B12) 



from which I find that can be written in the compact form: 



7g {w^ - K^a^) {w^ - K^v\) 



(B13) 



where w'^ = {lo — ku)^ , k? = {k — cou/c^)'^ , and where the fast magnetosonic speed perpendicular 
to the magentic field is given by (e.g., Vlahakis &; Konigl 2003) 

_ r 2 I 2 2 2 / 21 1/2 r 2 / 2 I 211/2 

It is easily seen that this expression for reduces to the relativistic pure fluid form 



2 f ''> ''2 



{ku — jj)^ 



,2\2 



{k — Loujc ) 



given in Hardee (2000) and that this expression for 0^ reduces to the non-relativistic MHD form 



P" 



n^a^) ( 



(a^ + v\) tu^ — K'^v\a'^ 



{ku — uj) 



(a2 + V^) {ku 



CO 



k'^VW 



k' 



where — > k'^ and v'j^ — > given in Hardee, Clarke & Rosen (1997). 

The solutions that are well behaved at jet center and at infinity are P*i{r < R) = CjJ±n{Pjr), 

and Pei{r > R) = CeH^l^{Per), respectively, where J±„ and are the Bessel and Hankel 

functions with arguments defined as 



P 



(B14a) 
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and 



7eM 



We 



1)^ Tip' 



(B14b) 



where = (a; - kuj^^f, = {k - uuj^e/c^f. lie = (l - ^iJc^) ^ and -f\^ = (l - v\J^^ \ 



The jet flow speed and external flow speed are positive if flow is in the +z direction. 

The condition that the total pressure be continuous across the jet boundary requires that 



The first derivative of the total pressure is given by 

d 



(B15) 



-K 



and with 



Vrl 



1 1 ■ 



-iXv, 



= —i{iO — ku) 



where is the fluid displacement in the radial direction it follows that 



dr 



{OJ — ku) X^r 



(B16) 



The radial displacement of the jet and external medium must be equal at the jet boundary, i.e., 
= (R) 5 from which it follows that 



f3j 



dJniPjr) 



-{u-kuj)Xj ^ d{/3jr) 



r=R 



— {u — kUe) X( 



-a 



(BIT) 



Inserting Cj and Cg in terms of the Bessel and Hankel functions leads to a dispersion relation 
describing the propagation of Fourier components which can be written in the following form: 



Pj j'niPjR) ^ fie H^n^' {peR) 
XjJniPjR) XeH^^\(3eR) 



(B18) 



where the primes denote derivatives of the Bessel and Hankel functions with respect to their argu- 
ments. The expressions 

Xj = l-lAjWj {zu] - k]v%) (B19a) 

and 

Xe = ihAeWe H ' l^lv^e) (B19b) 

readily reduce to the non-relativistic form x = Po[{'^ ~ ku)'^ — k'^Vj\ where Wq — > po given in 
Hardee, Clarke k, Rosen (1997). This dispersion relation describes the normal modes of a cylindrical 
jet where n = 0, 1, 2, 3, 4, etc. involve pinching, helical, elliptical, triangular, rectangular, etc. 
normal mode distortions of the jet, respectively. 
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C. Analytic Solutions and Approximations 

Each normal mode n contains a fundamental/ surface wave and multiple body wave solutions 
to the dispersion relation. The low- frequency limiting form for the fundamental/surface modes are 
obtained in the limit where u — > and k — > but with oj/k ^ In this limit the dispersion 
relation for the fundamental (n = 0) and surface (n > 0) modes is given by 



ra = 



(CI) 



and 



Xj 



-Xe n>Q 



(C2) 



where in this limit (ieR — ^ and fijR — > 0, and I have used the small argument forms for the 
Bessel and Hankel functions to write 



Um H^r^\p,R) 

where e is Euler's constant. 



-\{(5eR) {(3jR) [ln(%^) + f e - if 



n = 
n > 



C.l. Fundamental Pinch Mode (n = ; m = 0) in the low frequency limit 



In the low frequency limit, dispersion relation solutions for the fundamental axisymmetric 
pinch mode are obtained from equation (CI) 



r:' 



ln( 



(C3) 



Here we have the trivial solution tu| — n'^jv\j = with = v\j and the more interesting zeroth 
order solution 



msj 



with wave speed in the proper frame given by 



2 / 2 



31 J 



v%a] 



msj 



(C4) 



(C5) 



To first order this magnetosonic wave solution (eq. C3) can be written as 



vj] [1 - (5] ~ n]a] 



(C6) 
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where 



2 



ln( 



(C7) 



and S is complex. Thus, in the low frequency limit this fundamental pinch mode (n = 0) solution 
consists of a growing and damped wave pair with wave speed in the observer frame 



k 



Uj ± Vu, 

1 lb VyjUj/c^ 



(C8) 



where 



2 2 



1)2 



msj 



(C9) 



Previous work has shown the unstable growing solution associated with the backwards moving (in 
the jet fluid reference frame) wave. 



C.2. Surface Modes (n > ; m = 0) in the low^ frequency limit 

In the low frequency limit the fundamental dispersion relation solution for all higher order 
modes (n > 0) is most easily obtained from equation (C2) written in the form 



Te 



(CIO) 



where I have used x = IqIa^o — = 7o Wq [(a; — ku)'^ — (A;^ — u'^/c^) V|/7o] • The solu- 

tion can be put in the form 



UJ 

k 



['r]Uj + Ue] ± iri'^''^ {u^ - - VlJl]ll 



1/2 



(i + yl/7|c2) + r?(i + yi^./7f 



where 



and a "surface" Alfven speed is defined by 



_|_ 



ATrWjWe 



The jet is stable to higher order n > fundamental mode perturbations when 



7|7e' {n, - Uef < 7l,7le {Wj/lle + Weh%) ^^^^ 



Equation (Cll) reduces to the relativistic fluid expression 

UJ rjUj + Ue ir]^^"^ 



k 



1 + ri l + rj 



{Uj - Ue 



(Cll) 



(C12) 



(C13a) 
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given in Hardee & Hughes (2003) equation (6a) where for pressure balance and equal adiabatic 
index in jet and external medium r] — > jjae/jgaj. Similarly equation (Cll) reduces to the 
non-relativistic MHD form 



1/2 



1 + 77 1 + ?7 



(uj-Ue) -V; 



As 



1/2 



(C13b) 



given by Hardee & Rosen (2002) eq. (4) where V^^ — > {pj + pe) ^-B? + B^^ j {Ampjp^ and j] 

Pj/Pe- 



C.3. Body Modes (n > ; m > 1) in the low frequency limit 

In the low frequency limit the real part of the body wave solutions can be obtained in the limit 
a; = 0, A; 7^ where the dispersion relation can be written in the form 



cos [PjR - (2n + l)7r/4] 



_ XePe f 7rf3jR\ 



Xj Pj 



1/2 



^1 



(C14) 



Here I have assumed that the large argument form J„{3jK) ~ (2/7r/?ji?)^^^ cos [PjR — (2n + l)'7r/4] 
applies. In the absence of a magnetic field and a flow surrounding the jet, Xe = 0, e„ = 0, and 
solutions are found from PjR — (2n + l)7r/4 = ±m7r it 7r/2, where m is an integer. Provided 
e„ << 7r/2 and 6 k, cos^^ c„ ~ ± (7r/2 — e„), solutions can be found from ^jR — (2n + l)7r/4 = 
zb [mTT + (7r/2 ± e„)], where for ibe„ the plus or minus sign is for m odd or even, respectively. In 
the limit a; = 



(3^R 



1/2 



kR , 



(C15) 



and the solutions are given by 



kR: 



.mm o — 



7;2 o<2 



1/2 



X [(n + 2m - l/2)7r/2 + (-l)'"e„] 



(C16) 



where I have set m — ^ m + 1 to be consistent with previous notation so m = 1 corresponds to the 
first body mode. 

In the limit » equation (C16) reduces to the relativistic purely fluid form found in 
Hardee & Hughes (2003) 



^R ~ ^mn-^ = 



[(n + 2m - l/2)7r/2 + (-l)™e^ 



7j 



1/2 



(C17a) 



where M? = uj/aj. In the limit v\j » Oj equation (C16) becomes 

[{n + 2m - l/2)7r/2 + (-l)™e^ 



kR: 



imm Tj _ 



7j 



Ml, 



1/2 



(C17b) 



where M^j = uj/v^j. 
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Equation (C16) reduces to the non-relativistic MHD form found in Hardee &; Rosen (2002) 



1-MIJM%M] 



-1/2 



X [(n + 2m - l/2)7r/2 + (-l)'"e„] 



where = Uj/{aj + v^j) and I have used 

'j]{uj-a])iu^^-v%) 



1 - m2 ^/m1,m2 



1,2 . 



7.2 . 



M'i .M'4 

A] 3 



(C17c) 



We note here that there is an error in equations (5) in previous articles in the treatment of the sign 
on En for even values of m. 



C.4. The Resonance Condition 

The resonance conditions are found by evaluating the transmittance, T, and reflectance, 71, of 
waves at the jet boundary where T = 1 + TI. With the dispersion relation written as 



1 



1 Hi^'^'iPeR) 



where Z = x/P with 



(C18) 



.1/2 



the reflectance 



n 



(C19) 



(C20) 



iZe-Zj)/iZe + Zj) . 

For a fluid containing no magnetic field Z is a quantity related to the acoustic normal impedance 
(Gill 1965). When Zg + Zj 0, TZ and T are large, and the reflected and transmitted waves have 
an amplitude much larger than the incident wave. 



C.4-1- The Fluid Limit (Alfven speed <S sound speed) 
For the case of a pure fluid 



and 



Ce 



Cj 



(C21a) 



(C21b) 



where x/k'^ = W [C^ + j^K^/k'^] and C, = f3/k. For non-relativistic flows where {v? /c^){uj /ku) « 
1, 7 1, and with adiabatic indicies Vj = Fg the reflectance 



7^ = 



(Ce-0)(CeCi-l) 
(Ce + 0)(CeO + l) 



(C22) 
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and a supersonic resonance (Miles 1957) occurs when /3e + Pj = k{Ce + 0) = 0- This supersonic 
resonance corresponds to the maximum growth rate of the normal mode solutions to the dispersion 
relation. 

I now generalize the results in Hardee (2000) to include flow in the external medium relative 
to the source/observer frame. Here + Zj = becomes 



A necessary condition for resonance is (j < and Ce > 0; and on the real axis 



U-i — 



00 

> -r > 



Ue + Oe 



1 — Ujaj/c^ k 1 + Ueae/c^ 

It follows that the resonance only exists when 



u-i — a-i 



> 



1 — Ujaj/c^ 1 + Ueae/c^ 



or equivalently 



> 



(C23) 
(C24) 

(C25a) 
(C25b) 



To find the resonant solution for the real part of the phase velocity I solve C| = e^(^ where 



here I set e = (Tjj'jw^/k'^aj)/ {Vejgwj/k'^al) = 1 so that 



{uj/k — UjY 



1 



k C2 



{uj/k - UeY 



k 



(C26) 



The resulting quadratic equation can be written in the form 

v2 



where I have used 



+ 



C4 



(u]-afi-^la]{ul-afi 



(C27) 



7e«' 1 1 



and 7^ = (l — a^/c^) ^ . The solutions to equation (C27) are given by 







7se7sj J ^ 


— 2UjUe + "Ug 
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(C28) 



with the resonant solution given by 



"- = 1 = 



W {jseaehjUj + (-fsjajheU, 



(C29) 
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Inserting the resonant solution (eq. C29) into the expression for e gives 

0.695 <e^ = < 1.44 

else 

where 2.78 < r^7^ < 4. When Oj = Og and Tj = Fg, = 1, and the resonant solution is exact. The 
small range on e (0.83 < e < 1.2) suggests that this solution remains relatively robust for unequal 
values of the sound speed and adiabatic index in the jet and external medium. In the absence of 
an external flow the resonant solution 

/ \ 



■2^1/2 



k 7.-(7.eae) + ilsja,) _ V2 ^ ^^^^ _ ^ 

is equivalent to the form given in Hardee (2000). 

The resonant frequencies can be estimated using the large argument forms for the Bessel and 
Hankel functions. In this limit the dispersion relation becomes 

X(ftfi)£(Afi) ^ _ 2n + l ^ X, A (C30) 

From the dispersion relation with Zg + Zj ~ 0, and {Xj/Pj)iPe/Xe) = ^j/^e ~ Pe/Pj ~ — 1, 
tan[f3jR — (2n + l)7r/4]Re ~ on the real axis. It follows that 

\l3jR\ « \peR\ ~ (2n + l)7r/4 + rmr 

can be used to obtain an estimate for the resonant frequencies from |/3e-R| ~ (2n + l)7r/4 + rmr, 
with result that the resonant frequencies are given by 

^nmi? _ l^nmR ^ (2n + l)7r/4 + 17177 ^^^^^ 

In the absence of external flow, Ue = 0, and for Uj » and 1 >> (/cag/w)^ this expression 
rediiccs to the form given in Hardee (2000). 

When 7j(7sgae) >> ^ei'JsjO-j), the resonant wave speed becomes ^ Uj, ng/v^ Ue/uj and 
provided ttg << uj and Og << uj, the resonant frequency increases with increasing u^/uj and a^/uj 
as 

^ {2n + l)7r/4 + mvr ^^^^^ 



l-2V^,(l-a2/c2)-(a2-u2)/u2 



1/2 



^J\^ "-e/ J \^e "-e)! "-j 

In general, the resonant frequency o;*^ — ^ oo as (1 — u^jv*^ — [a^jv^ — Ueae/c^)^ — ^ 0. An 
equivalent condition for (1 — Uf^jv^ — {a.e/v'^ — u^aej^^ = is 

tXj — Ug aj + Og 



l — UjUe/c^ l + ajOg/c^ 



(C33) 



and the resonance moves to higher frequencies with w*^ — > oo when the "shear" speed drops below 
a "surface" sound speed. 
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The behavior of the growth rate at resonance also can be found using the large argument forms 
for the Bcsscl and Hankcl functions. In this limit the reflectance can be written as 

K = Ml = MmHi!yM-j:(mHl!\m ^ - 5i±i.)i , (C34) 



and 



^ „ 2n + 1 i , ,„ , (t> 
/3,•i^-^^7r«.-ln|7^|-| 



where 7^ = |7^| e'"^. It follows that 



iPjR)i -ln\n\ 



and since typically at resonance, \u! — kRUj \ /aj > |A;_r — (jJUj/c^\ I can approximate Pj by 



It follows that 



and 



{iu - kRUj 



— iki 



1 a. 



■ln|7^| . 



At resonance 



7^ 



2 7jUj 

(Ze - Z,) _ - /3e(Xj/Xe) _ /?, " /?e 



-2/3e 



and 



1/2 



(C35) 
(C36) 

(C37) 

(C38) 

(C39) 

(C40) 
(C41) 



where I have used — f 
that 



2/3^ from the resonance condition on the real axis. It follows 



\n\ 



lj--le — 
'J ai ae 



-^^kR 



kj 



where I have used 



■(w - kRUe) 



1/2 



(C42) 



iki 



+ 



-kR 



then 



^Oe U - kRUe 

If I assume that 7j(7seOe) >> leilsjO-j), with resonant wave speed uj/kR ^ Uj and Ue/uj « 1 



\n\ 



4 (^^)^ (1 - 2Ue/Uj) + 



(1 + Me/lij) 



1/2 



(C43) 
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and since kjR ^ — {aj/2jjUj) In \TZ\ 



mi 



1/2 



(C44) 



From equation (C32) 



(1 — 2Ue/Uj 



(1 - 2Ue/Uj) 



l-2{ue/uj){l-al/c^)-al/u] 



[(2n + l)7r/4 + m7r] 



and if say Ue = 0, then 

(l^l' - l) 



,9 \V2 , , 

7^P-1| ln7^«4 



l-a2/u2 



1/2 



[(2ra + l)7r/4 + tott] 



(C45) 



Formally \TZ\ — > oo as a;*^ — oo when the jet speed drops below the "surface" sound speed 
given by equation (C33). This result applies only to the surface modes and not to the body modes 
as, in the fluid limit, the body modes do not exist when the jet speed drops below the jet sound 
speed, see equation (C17). On the other hand, if say, a^/uj « 1, then 



(^\nf - l) In \n\ « 4 [{2n + l)7r/4 + rmr] 



(C46) 



Formally \TZ\ « constant as o;*^ — > oo when the wind speed becomes comparable to the jet speed, 
Ue < Uj, as must be the case for the velocity shear driven Kelvin-Helmholtz instability. 



C.4-2. The Magnetic Limit (Alfven speed S> sound speed) 
For the magnetic limit in which magnetic pressure dominates gas pressure 



and 



(C47a) 
(C47b) 



Zj = I'j'yljWjVAj {-cu] - K^jV^f^ , 

A necessary condition for resonance is (tUg — Kg'^^g) > and — i^^v^j^ < on the real axis 
with result that Z^. + Zj = Q when 



Uj - VAj 



UJ Ue + VAe 

> - > 



(C48) 



1 — UjVAj/c^ k 1 + UgVAe/c^ 

It follows that the resonance only exists when 

Uj - Ue ^ VAj + VAe 
1 - UjUe/c^ 1 + VAjVAe/c^ 

This result is identical in form to the sonic case with sound speeds replaced by Alfven wave speeds. 



(C49) 
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The resonant solution for the real part of the phase velocity is obtained from 

(C50) 

where I have used 7^7^ [tu'^ — i^'^v\) = 7^ [zu'^ — [k'^ — co'^/c^) V\h'^\ , and recall that v\ = V\/^\. 
The resulting quadratic equation can be written in the form 



+ 



(C51) 



where I have used 



because pressure balance in the magnetically dominated case requires WjV^j = 
solutions are given by 

k ^jwjvl^ - ilwlvl^ 

and the resonant solution becomes 

* _ ^ _ li^jVAjUj + -jeWeVAeUe _ (7Aef Ae) ljUj + {nAjVA:]) le^e 



W,Vl^. The 



k 



Ij ilAeVAe) + 7e {lAjVAj) 



(C52) 



(C53) 



where I have used WVa = WVj^/ {'yAVA)j and WjV^- = WpV^^. This resonant solution has the 
same form as the sonic case with sound speeds and sonic Lorentz factors replaced by Alfven wave 
speeds and Alfvenic Lorentz factors. 

As in the sonic case the resonant frequencies are found from 

\PeR\ ~ (2n + l)7r/4 + mn 

with result that the resonant frequencies are given by 



(2??, + l)7r/4 + m7r 



VAe 



VAe 



(1 - Ue/vl,f - {vAe/v^ - UeVAe/c^f 



1/2 



When 7j {'jAeVAe) » 7e ilAjVAj) the resonant wave speed becomes ^ uj and Ue/vl 
and provided Ug « Uj and VAe « Uj the resonant frequency increases with increasing Ug/uj and 
VAe/uj as 

(2n + l)7r/4 + m7r 



VAe 



1 - 2 {uju^) (1 - v\Jc^) - {v\^ - ul)/u] 



1/2 • 



(C54) 

nd 
(C55) 



Here the resonant frequency w*^ — > 00 as {1 — Ue/v^j)'^ — [vAe/v^ — UeVAe/c^Y 
equivalent condition for (1 — Ue/vljf' — {vAe/vl^ — u^VAe/d^)'^ = is 

Uj -Ue _ VAj + VAe 



1 - UjUe/c^ 1 + VAjVAe/c 



,2 ' 



0. An 



(C56) 
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and the resonance moves to higher frequencies as the "shear" speed becomes trans- Alfvenic. 

The behavior of the growth rate at resonance proceeds in the same manner as for the fluid 
hmit but with sound speeds replaced by Alfven wave speeds. The resonant growth rate is now 
given by 

1 VAj 



kjR. 



2 7j-u, 



In |7^| 



(C57) 



If I assume that ^j{'~^AeVAe) » leilAjVAj)-, with resonant wave speed uj/Ur^ uj and Ue/uj « 1 
then 



VAe 



(1 - 2ue/uj) + kjR^ 



7,^-^(l + «e/n,) 



Prom equation (C54) 



[1 - 2Ue/Uj) 



(1 - 2Ue/uj) 



1/2 



(C58) 



VAe 

and if say Ue = 0, then 



1 - 2 {Ue/Uj) (1 - VIJC^) - {V\^ - Ul)/U] 



[(2n + l)7r/4 + rnvr]' 



{\n\^ - 1) 



1/2. 



ln|7^| 



1 - -ie/^j 



1/2 



[(2n + l)7r/4 + mvr] 



(C59) 



and \R.\ increases as a;*j„ increases when the jet speed decreases. However, when the shear speed 
drops below the "surface" Alfven speed, see equations (Cll & C12) the jet is stable. This result 
is quite different from the fluid limit where the jet remains unstable when the shear speed drops 
below the "surface" sound speed. If I insert 



1 — UjUf.ji? 
1 + VAjVAe/c^ 



{vAj + VAe) ■ 



from equation (C56) into equation (C12), it follows that the jet will be unstable when resonance 
disappears only when 



2 2 

.,2.,2 ^1 „. /„2^2 ^ o.,2 .,2 ^Ae + ^Aj ,„..„.. /„2\2 



2 2 /i / ^ n 2 2 



, , ,2(1+ VAjVAe/c ) . 
[VAj + VAe) 



(C60) 



where I have used [vl^ + v%) = [Wjh\^ + VFe/7i, ) \^] + ) / (4vrI^jVF, 



as Bf. = Bj from 



magnetic pressure balance. Formally \R\ — > 00 as w*^ — > 00 only for jet Lorentz factors greatly 
in excess of the Alfvenic Lorentz factor. 



C.5. Wave modes at high frequency 

To obtain the behavior of wave modes at high frequency I begin with the dispersion relation 
written in the form 
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and assume a large argument in the Hankel function with Hn\PeR) ~ expz [PeR — (2n + 1) 7r/4] 
and a small argument (5jR « 1 in the Bessel function to write 



PjRXe 



Jo(fijR}_ in/2 „ _ n 

Me-W2 „ > 



The small arguement form for the Bessel function gives Jo{(3jR)/Ji{PjR) ■^2/f3jR with result that 
the dispersion relation becomes 



Xj 



1/2 



-in/2 



n = 
n > 



(C62) 



At high frequency and large wavenumber Xj and Xe, are proportional to k'^, Pj and Pe are pro- 
portional to k, and PjR = CjkR oc {kR)^^"^ for n = 0. Thus, the internal solutions in the high 
frequency and large wavenumber limit are given by ~ and are found from 



{kuj — (j)^ — {k — LOUjl(?^ Oj {kuj — u})^ — {k — ijJUj/(?Y Vaj 
and it follows that 



„2\2 „,2 



UJ 



Uj ± aj 



or 



CO 



k 1 ± Ujaj/c^ ' 

Uj ± VAj 



k l^UjVAj/c^ 



(C63) 
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(C64b) 
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